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Abstract 

We study the structure of the manifold of sohtary waves in a particular three-component scalar 
field theoretical model in two-dimensional Minkowski space. These solitary waves involve one, two, 
three, four, six or seven lumps of energy. 

1 Introduction 

The search for solitary waves is a mandatory topic in a huge number of branches of non-linear science. 
The reason for this lies in the fact that the method of linearizing nonlinear differential equations omits 
phenomena of essential importance in many different physical problems. The non-linear Klein-Gordon 
system of coupled PDE equations, generalized to n-real scalar fields, reads: 

#-^^*" + S:=° <.=i.2....,« . (1) 

Here, l^((/>i, . . . , 4>n) is a differentiable potential function of the scalar fields (j)a, whereas are non-linear 
functions of the fields. It is clear that this system is invariant under Poincare transformations. For this 
reason we shall use a system of units where the speed of light is set to one: c = 1. Also, due to the 
potential applications of the PDE system in quantum physics we shall implicitly understand h = 1, 
i.e., we shall work in the natural system of units. 

We shall focus on physical systems governed by the PDE system in only one spatial dimension. 
The paradigm with only one real scalar field is the celebrated sine-Gordon equation: 

^ ~ 7^ + sm = 

Arising in differential geometry in the theory of negative curvature surfaces, the sine-Gordon equation 
is ubiquitous in the description of one-dimensional physical phenomena. To mention a few, 1) it is the 
evolution equation for the amplitude of various slowly varying waves, 2) it describes the propagation of 
a dislocation in a crystal, 3) it provides a model for elementary particles moving on a line, 4) it rules 
the propagation of magnetic flux in a long Josephson-junction transmission line, 5) it determines the 
modulation of a weakly unstable baroclinic wave packet in a two-layer fluid, etcetera, see, e.g.. Reference 

Point 3) deserves special mention and will be the arena within we shall analyze the generalized non- 
linear Klein-Gordon PDE system proposed in this paper. The similarities between the properties of 
solitary waves (or solitons) and elementary particles are clear. The field energy density of a non-linear 
wave (solitary wave or soliton) is localized in a lump. Moreover, being non-dispersive, solitary waves 
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(or solitons) propagate without changes in shape. After colhsions, sohtary waves almost keep their 
structure, leaving some radiation in the form of dispersive waves as remnants. Solitons, however, only 
suffer phase shifts under collisions. In soliton-antisoliton scattering, either both solitary waves may be 
annihilated or, alternatively, they may form a bound (breather) state. Elementary particles share all of 
these properties. Thus, if an appropriate system of non-linear field equations admits solitary waves (or 
solitons) as solutions, these non-dispersive waves behave almost as elementary particles. 

Nevertheless, the main feature of these solutions, the confinement of the energy density, has straight- 
forward implementation in many areas. In Condensed Matter, these solutions describe interfaces in 
magnetic materials j2j and in ferroelectric crystals [H]. They have also been used to understand some 
bizarre properties of the poly(oxyethylene) |^, widely studied owing to its biotechnical and biomedical 
applications 0. Among the impressive features of this polymer, we can mention the ability to entrap 
metallic ions in aqueous solution and its solubility in water to almost any extent. On the other hand, 
in Cosmology solitary waves are interpreted as domain walls which seeded the formation of structures 
in the early universe They are also of direct interest in high energy physics, for instance, describing 
gravity in warped spacetimes involving D extra spatial dimensions [7j, or the coupling of scalar matter 
fields with dilaton gravity [H], which seems to be related to the formation and evaporation of black holes 
[5]. Furthermore, when these defects appear as BPS states in both extended supersymmetric gauge 
theories |lflj and string/M theory they play a crucial role in the understanding of dualities between 
the different regimes of the system. In this framework, they behave as extended states in = 1 SUSY 
gluodynamics and the Wess-Zumino model [12]. The structure of the solitary waves that we are going to 
unveil in this paper is richer than the usual kink structure found in simpler models. Thus, our topological 
defects could describe more subtle effects in each of all these scenarios. 

In general, the above mentioned systems involve a high number of fields. In order to investigate the 
presence of solitary waves or topological defects, the usual procedure is to obtain an effective scalar field 
theory, carrying out severe restrictions on the original theory. In most cases one is compelled to pursue an 
effective theory that corresponds to a single scalar field model, where the existence of topological defects 
can be checked easily. In general, however, the effective theory depends on several scalar fields, and the 
truncation can involve a important loss of information about the presence of solitary waves or topological 
defects. Therefore, it is desirable to investigate the general properties of solitary waves in a multi-scalar 
field theory. This is an important qualitative step as reported by Rajaraman This already brings 

us to the stage where no general methods are available for obtaining all localized static solutions, given 
the field equations. However, some solutions, but by no means all, can be obtained for a class of such 
Lagrangians using a little trial and error. 

Some work on two-component scalar field theory models equipped with a Minkowskian space-time has 
been accomplished, see for instance References ^1 El El El El El EOl 1^ • In some of these models, 
two-parametric families of solitary waves or kinks are identified. Generally, each member of these families 
is a composite solitary wave, which involves several lumps, i.e, the energy density is localized at several 
points. The parameters identifying each solution correspond to a translational parameter, xq, which 
determines the center of the kink, and an "orbit" parameter b, which specifies the separation between 
the lumps. Study of Manton's adiabatic motion of solitary waves |22) is very simple in this case because 
the geometric metric of the moduli (parameter) space that governs the dynamics does not depend on 
the translational parameter xq; thus, it is always possible to apply a transformation that leads to an 
Euclidean metric. On the other hand, there are only a few works that have addressed research into 
solitary waves in three-component scalar field theory models, see |231 124| l25| I26j . In this richer case there 
are three-parametric families of kinks, and the adiabatic evolution of three-body lumps is associated with 
metrics with curvature. The low-energy dynamics of, in this case, three-body solitary waves is therefore 
much more intricate. 

In a slightly different physical context, where the solitary waves are understood as domain walls 
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(2-branes) in (3+1) dimensions models with three scalar fields have been discussed in Reference |27| . 
Considering topological defects (0-branes) in more spatial dimensions (p-branes) allows for complex 
structures having defects nested inside defects. It is plausible that the composite solitary waves to be 
discussed in this paper will give rise to even richer structures of nested defects inside defects if considered 
in more dimensions. An interesting work is also available on networks of topological defects in a model 
with two complex scalar fields, see The potential energy density is a polynomial in the fields of 

arbitrary order chosen in such a way that the model enjoys a U{1) x Z„ symmetry, networks of domain 
walls arising in the surface of Q-balls. In other direction, defects inside defects in models with two scalar 
fields have been considered before, see (22] and jSOl- Very recently, networks of topological defects arising 
in a model with 9n scalar fields in nine spatial dimensions have been studied in [HJI and the evolution of 
these (or similar) structures has been analyzed in j32j . 

In this paper, we shall address a three-component scalar field model that generalizes the two- 
component model discussed in Reference |14j to three fields. The organization of the paper is as follows: 
in Section 2 we introduce the model and describe the spontaneous symmetry breaking pattern as well as 
the spectrum of dispersive waves. In Section 3 we describe the variety of solitary (non-dispersive) waves 
that arises in this system. We classify these solutions according to the number of lumps they are made 
of. There are three basic solitary waves with the energy density confined in one finite region. The rest of 
the solutions are composite and display several lumps -the energy density is localized at various points-. 



2 Deformation of the A(0 • 0)2-model 

We shall deal with a three-component scalar field theory model defined in a (1+1) Minkowskian space- 
time. The dynamics is governed by the action functional: 



S= d^x 



(2) 



and the potential term is the following sixth-degree polynomial in the fields: 

v{$) = {cfi + 4>l + (t>l){(t>l + 4 + 4- 1)' + '^(44 + 44)i4 + 4 + 4-'^) + 44 + 44 ■ (3) 

£72 and (73 are the non-dimensional coupling constants of the system, chosen, without loss of generality, 
such that (72 < CTs . 0" , a = 1,2,3, are dimensionless scalar fields, 



4>{xf),xi) = {(i)i{xQ,xi),(j)2{xf),xi),(j)z{xf),xi)) : 



Our convention for the metric tensor components in Minkowski space M^'^ is: goo = ~gii = 1 smd 
gi2 = 921 = 0. 

This model is a (ct2, (T3) deformation of the same system with 50(3)-invariant potential energy density: 

vs{$) = {4 + 4 + 4){4 + 4 + 4- 1)' (4) 

Besides the Poincare transformations acting on E}'^, the non-deformed model -(T2 = cr2 = 0- is invariant 
under 0(3) rotations in internal space M^. The specific form of Vs, however, shows spontaneous symmetry 
breaking of the internal symmetry to a SO (2) subgroup so that two Goldstone bosons are unavoidable 
if the asymmetric vacuum -any point at a distance 1 from the origin in M^- is chosen to perform the 
quantization of the system; the set of zeroes of Vs is the the union of a discrete point, the origin, and a 
continuous manifold, the 2-sphere of radius 1 in M^. Coleman proved in [HH] that there are no Goldstone 
bosons on the line in a sensible physical system. The infrared asymptotic behavior of the quantum system 
would require modification of the potential Vs in such a way that the zeroes of the potential become a 
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discrete set and massless particles are forbidden. Our deformation complies with this requirement and 
could be understood as being generated by infrared quantum fluctuations. The added terms in the 
potential © explicitly break the 0(3) symmetry to the discrete subgroup G = Z2 x Z2 x Z2, generated 
by the reflections (pi — > —(pi, (p2 —<p2 and (p^ — > —(p^ in internal space. We also stress that the present 
system is a generalization of the two-component scalar field theory model discussed in Reference |14) . 
where the dimension of the internal space, the number of scalar fields, is increased by one unit. 
The field equations form the system of coupled second-order partial differential equations: 

^'"^^ - 2(Pi[3icpl + 4 + 4f + ^-^4-m-4)4-2{2-4y'^ 







d^(p2 






dx'^ 


d'(Ps 






dx"^ 



^31 



2(P2 [3(</.? + (Pl + (Plf - 2(2 - al)(Pi - \al(g - 2{al + (jI)cPI + 4] (5) 
203 [3(0? + (Pl + <Plf - 2(2 - al)cPl - 2{al + al)(Pl - AalcPl + 4] 



For the sake of simplicity, henceforth we use the notation: = t, xi = x, 02 = yl^-frf, <73 = y^l — o". 

(Ton = a\- (To. 



2.1 Structure of the configuration space 

A "point" in the configuration space of the system is a configuration of the field of finite energy; i.e., a 
picture of the field at a fixed time such that the energy E, the integral over the line of the energy density, 

1 f d(Pi \ 2 1 / d(p2 \ ^ 1 / d(p'i ^ 



is finite. Thus, the configuration space is the set of continuous maps from M to M'^ of finite energy: 

C = |0(x) G Maps(M,M^)/£; < ooj 
In order to belong to C, each configuration complies with the asymptotic conditions 

lim (PeM lim — = (7) 

where is the set of zeroes (minima) of the potential term V{(p). In our model, Ai consists of seven 
elements, see Figure 1: 



M 



0° = (0,0,0); 0^* = (±1,0,0); 0^* = (0, 0, icxg); 0^* = (0, ±^2, 0)} (8) 



Note that all the zeroes except the origin lie in the ellipsoid ]E = 0i + ^ + ^ = 1- The intersections 



^2 ^2 



of the ellipsoid E and the 0i = 0, 02 = 0, 03 = 0-planes are respectively the ellipses ei = ^ + 

2 2 2 ""2 "3 

62 = 0f + % = 1 and 63 = 0f + % = 1 of eccentricities e{ei) = £(62) = o"? and £(63) = 4- The 
minima B± are on the 03-axis at the intersection between ei and 62; C± lies on the 02-axis at ei H 63, 
and D± on the 0i-axis at 62 H 63. 

The configuration space C is the union of 49 topologically disconnected sectors: 

c = ui,jC'^ 

where I, J = {O, D, C, B} and C^"^ stands for the set of configurations that connect the points (p^ with 
0"^. Because temporal evolution is a homotopy transformation, the asymptotic conditions do not 
change with t and the sectors are completely disconnected. Physically, this means that a configuration 
in a sector cannot evolve into configurations in other different sectors; it would cost infinite energy. 
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A3 

Figure 1: Points belonging to Ai, umbilical points, focal points, and characteristic curves of the model drawn on 
and inside the ellipsoid ¥,: (left) Cartesian coordinates, (right) Elliptic coordinates. 

2.2 Plane waves and spontaneous symmetry breakdown 

The zeroes of the potential energy density are also homogeneous static solutions of the field equations 
^ with zero energy; in fact, they are absolute minima of the energy in the topological sectors C^^ . A 
spontaneous symmetry breaking arises because of the degeneracy of the set A4. Whereas the point (j)^ 
located in the origin preserves the discrete symmetry Z2 x Z2 x Z2 of the Lagrangian, the choice of cj)^ 
breaks the symmetry generated by the reflection cfti —(t>i, which connects the two points, (j)^^ and 
(j)^ , to the little group {e} x Z2 x Z2. Similar considerations work for the pair of points and cj)^ 
with respect to the transformations (j)2 —(1)2 and (p^ —(ps- The moduli space of vacua comprises four 
G-orbits 

where we write 4>^ = {(j)^^ , (j)^' } , = {(f)^^ ,4f~} and cj)^ = {(^^^,(^^~}. 

Small fluctuations around the homogenous static solutions (f), ^^{t, x) = (p + Sip^t, x) are still solutions 
of © if the linear PDE system 

{D6a, + m2,(<^)) 6Mt, x) = Oa, Ml = ^£^(0) (9) 

is satisfied. The solution of © via the separation of variables 5ipa{t,x) = Re (e*'^*) f'^{x) leads to the 
spectral problem for the second-order fluctuation -or Hessian- operator: 



6=1 6=1 ^ ^ 



x)=a;^/-(x) (10) 

H(0) is a diagonal matrix differential operator for the four points of M with a positive definite spectrum; 
each constant solution belonging to is stable. Thus, there are four types of dispersive wave-packet 
solutions, each one with three branches 

/r(x)=sinA:x,/2" = /3- = ; /- = Q , /^-(x) = singx , /g^ = ; = = , (x) = sinpx , 

living in one of these sectors. The building blocks are these plane waves and the dispersion laws of the 
corresponding wave packets are respectively determined by the diagonal matrices: 





2a\ 
M'^{^^') = ( 2f7f2 
8a| 
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The three branches for each type are: 













Cj2(A;) = A;2 + 2(j| 

i^Hq) = q' + 8al 

io\p)=p^ + 2al^ 


a;2(g) = g2 + 2al^ 
coHp) =P^ + 8al 


u;^{k) =e + 2 

u:\q)=q^ + 2al 
^Hp) =P^ + 2al 



2.3 Solitary waves from integrable dynamical systems 



According to Rajaraman [T3] "A sohtary wave is a locahzed non-singular solution of any non-linear 
field equation whose energy density, as well as being localized, has space-time dependence of the form: 
£{t,x) = £{x — vt), where v is some velocity vector". Because we are dealing with a Lorentz invariant 
system, solutions with a temporal dependence like this are obtained from static solutions by means of a 
Lorentz velocity transformation: L{\)(j){x) = 4>^{t,x) = ^ ' 



Thus, the search for solitary waves is tantamount to looking for localized solutions in the ^ J 

sector of the ODE system : 

Understanding the scalar field (p ^ the particle position vector; x as the particle time and U = —V 
as the particle potential, this is no more Newton equations describing the motion of a particle in three 
dimensions in the field of the potential U. Moreover, the energy functional and the energy density © 
of the field theory turns respectively into the action and the Lagrangian of the mechanical system. The 
Legendre transformation provides the Hamiltonian of the analogous mechanical system: 

a=l a=l 

The Hamiltonian l\12\\ corresponds to a Stackel separable system. To prove this, we introduce three- 
dimensional Jacobi elliptic coordinates {X,fi,u) in the internal space {(pi, (1)2, (ps), see [22]. In terms of 
elliptic coordinates, the fields are: 

) 

(13) 



(l-^32)(l-a2) ' ^2 (^2_^2)(^2_i) ' ^3 (^2 _ ^2) (^2 _ i) 



The new variables take values in the infinite cuboid 

— OO <A<fT|<^<fT|<l^<l 

which is mapped by the coordinate transformation ()13|) to a Cartesian octant in M^. Therefore eight 
points in become only one in the cuboid and mapping back solutions in elliptic coordinates to the 
original Cartesian ones requires a lot of care to reassemble the solution in the total internal space from 
one octant by imposing continuity of the trajectories and their derivatives. Apparent rebounds on the 
faces of P{0), see Figure 1, are really continuous in the closed domain in bounded by the ellipsoid E. 

In Figure 1 we have plotted the finite elliptic cuboid, denoted by -P(O), where (an octant of) the 
interior of the ellipsoid E is mapped. The ellipsoid itself is mapped to the A = face of P{0), E = A = 0. 
In elliptic coordinates, the following information is interesting: 

• Homogeneous solutions 

4>^ = {X = 0,^ = aliy = l) , = {X = 0,fi = alu = l) 
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Ellipses where the homogeneous solutions live: 



A = , = cj| 
A = , V = al 



ei = A = 0,i/=l ; e2 = <{7_n'!^_l2 ! = \ = Q , ix = al 



• Other special points: Fi, F2, -F3 (foci of ei, 62 and 63 respectively) and the umbilicus points A of 
the ellipsoid E are settled on the vertices of the parallelepiped P(0) in elliptic space: 

Fi = {\ = al,^ = al,v=l) , F2 = {\ = al,n = al,v = al) 

F3 = (A = = a|,z^ = , A= {\ = 11 = al,v = al) 

4,2 j,2 



Special curves such as the ellipse 64 = ^ H — |- = 1 corresponding to the edge F1F3 in P(0) and 

22 "^^ "^^^ 

the hyperbola /i = % — ^ = 1, the edge F2A. 



'2 



ei= \ = al,v = a\ , h= y, = a\,v = a\ 

We shall see that these curves are focal lines of finite action trajectories in the analogous mechanical 
system. 

In elliptic coordinates, the Hamiltonian (|12() of the analogous mechanical system becomes: 

^ ^ -^A ^ ^ ^ ^^^^ 

{X-fi){X-u) {fi-X){fi-u) {u-X){u-fi) 

Here, 

Hx = 2P3 {X)pI + P4 ( A) ; H^ = 2P3 {fi)pl + P4 (p) ; = 2P3 {u)pl + P4 {v) 

P^{x) = (1 - x){al - x){al - x) Pi{x) = x^{x - l){x - al){x - af) 

Expression (|14)) is of the Stackel form. Thus, the analogous mechanical system is Hamilton-Jacobi 
separable. We will take advantage of this fact in the following Sections in order to identify the variety of 
solitary waves in the model as the finite action trajectories of the analogous mechanical system. 



3 The variety of solitary waves 

In this Section we describe the variety of solitary waves in successive stages. The analogous mechanical 
system encompasses a vast manifold of finite action trajectories in one-to-one correspondence with the 
solitary waves of the field theoretical model. One-body solitary waves, made of only one lump, can be 
found by means of Rajaraman's trial orbit method, without the need to use elliptic coordinates. On 
plugging some specific trial curves into the equations simple solitary waves that we will refer to as 
basic lumps arise. The trial orbits correspond to some edges of the parallelepiped -P(O) shown in Figure 
1(b). Later, solutions which comprise several of the basic lumps will be found by restricting the movement 
to the faces of the parallelepiped -P(O) to find solitary waves akin to those discovered in Reference ^1] 
. As a last step, we shall solve the equations ((TT|) in the general case by applying the Hamilton-Jacobi 
theory with no restrictions. 

The system of ODE (jlip is invariant under translations x ^ x — xq and reflections x —x in the 
spatial parameter x. The first symmetry allows us to localize the center of the solitary wave or kink at 
an arbitrary point xq, whereas the second symmetry connects a kink with its antikink. Henceforth, we 
shorthand these transformations on x as x, i.e, x = (-l)'^(x - xo) with xo G M and (5 = 0, 1. 
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3.1 Solitary waves: basic types 



1. K'^^ : The solutions of (fTT|) in the sector C^^ , located in the (/>3-axis in internal space, are four-fold, 
see Figure 2(c): 

02(x) = O; 03(x)= (Zl)!^ (15) 



(/)i(x) =0: 



Here, a = 0, 1 distinguishes between solutions with cj)^ > {) and (/)3 < 0. We refer to these solutions as 

ROB. 

Some remarks about this notation: the superscripts specify the elements of M that are connected 
by the solution and/or the topological sector C*^^, where the solitary wave lives. The subscript shows 
the number of basic particles, lumps or bodies of which a particular solitary wave is made. In Figure 
2(a) a a = (5 = K^^ solitary wave of kink shape is depicted; this K'^^^ kink connects the points 
-B_|_ and O m. }A and lives in the C^^^ topological sector of the configuration space. In Figure 2(b) 
its energy density is plotted as a function of x; because it is localized at a single point we understand 
these one-body solutions as basic (non-composite) particles. A Lorentz transformation applied to 
will show this lump of energy to be a solitary (traveling) wave. Integration on the real line provides us 

-4 

with the total energy, which for these solutions is E{K^^) = and this depends only on the coupling 
constant a^. 




0.1 
0.05 



-30 -20 -10 




Figure 2: Kink form factor (a), energy density (h) and orbits (c) for K^^ kinks. 



We shall adopt the policy of showing all the energy density plots with the same vertical range in order 
to compare the size of the lumps associated with different solitary waves. 

2. K^^ : Plugging the trial elliptic orbit 63 = 0f -|- ^ = 1 into the equations (fTT|) we find the eight 

""2 

solutions 



01 (x) 



02 (x) 



03 (X) = 



(16) 



Vl + e^v^^ Vl + e^ 

with a, 5 = 0, 1. These trajectories are embedded in the 03 = plane and live in one of the eight C^^ 
topological sectors, see Figure 3(c). Note that both non-null components, 0i and 02, have kink shape. 

The value of 5 determines whether the solution is a kink or antikink. The values of a and /9, however, 
characterize the quarter-ellipse in the intersection of the ellipsoid E with the 03 = plane, where a 
particular solution is located; therefore, a and (3 also choose the points that are asymptotically connected 
by these kinks. In Figure 3(a) a particular kink K^^ with a = [3 = 5 = Q \s shown that connects the C+ 
and Dj^ points. The energy of these kinks is E{K^'^) = "^^/^^j ^ function of the coupling constant o"2 
only. Again, the energy density is localized at one point, and they are basic particles in our model, see 
Figure 3(b). 

3. KP'-^ : We shall also try the elliptic orbit ei = % + % = 1 confined to the plane 0i = 0. 
Substituting this condition in (fTT|) . the following eight solutions, connecting B and C, are obtained: 



01 (x) =0: 



02 (x) 



(-1)"^2 



{-If 



(17) 
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Figure 3: Kink form factor (a), energy density (h) and orbits (c) for K^^ kinks. 



with a, (3,6 = 0, 1. Again the non-null components, 02 and (p^ in this case, are kink shaped. In Figure 
4(a) a Kf'-" kink is represented for a = (3 = 6 = 0. The analogous trajectory leaves the instability point 
B-^- U = —V "at" X = — oo and arrives at another instability point C+ "at" x = oo. These solutions 

-4 -4 

are made of a single lump and carry energy equal to E{Kf^) = '^^^J^ ' ^ function of both coupling 
constants. 




Figure 4: Kink form factor (a), energy density (h) and orbits (c) for K^'~^ kinks. 



In sum, there exist three kinds of one-body solitary waves or basic extended particles, denoted as 
and i^f^. The energy increases as the eccentricity of the ellipse 63 grows, whereas the 

energy of K'^^ decreases when the eccentricity of 62 increases. The energy for the K^'-' lumps depends 
on the difference of the eccentricity of the ellipses 63 and 62 and becomes zero when these magnitudes 
are equal. Which basic particle is more energetic depends on the values of the coupling constants 02 and 
(T3. The figures in this paper have been drawn for the values (T2 = 0.5 and 0-3 = 0.8. 



3.2 Two-body solitary waves 

1. {b): We shall now tackle the search for solutions in the (/>! = plane, v = I'm. elliptic coordinates, 
see Figure 5(c). The Hamiltonian (|14|) restricted to this plane reads 

2P;{\)pI + P^{\) 2P3(/i)p^ + A(M) 
(A-/i)(A-l) + (^-A)(/i-l) 

The standard Hamilton-Jacobi procedure prescribes how to solve the HJ equation 

dS [dS dS , \ 

assuming the separable form of the Hamilton's principal function S = Sx + S\{\) + 5^(;u). The PDE HJ 
equation becomes a ODE system and the complete integral can be found by quadratures 



-Ex\ Sx = Sx / d\ 



F-\EX 
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where 72 is a 



in terms of the integration constants E, F. Sj = Signpj denote the sign of the momenta. 

The trajectories are determined from the orbits selected by the equation 72 
constant, and the time-schedule fixed by the equation 7i = where 71 is another arbitrary constant. 
The asymptotic conditions picking finite action trajectories force E = F = 0. The finite action 
trajectories of the analogous mechanical problem are thus determined by the pair of equations 



,^272 



\-4 


-A" 










[-4 


-A. 











which can be reshuffled in the form: 

-^/2a2(2i-<7|72) 



A 


sa 











o-|72) ^ 


" A " 


Sa 






[^I-aJ 







(18) 



The choice sx = s^ in (|18|) yields solutions in the topological sector C*^*^. Translating (|18|) back to 
Cartesian coordinates, we obtain the following field profiles for the associated one-parametric family of 
solitary waves parametrized by 6^ = 



(-l)°a2 



1 _|_ "32 ^-2y/2a^x _|_ Z2 ^2g-2V2(T|2S 
"3 



52 _|_ '^32g-2V2CT|x _|_ ^3 g2v^(T|2S 



. (19) 



We refer to these solutions as K2^{b) because they are made of two basic lumps, a = describes 
solutions joining O and C'^ whereas a = 1 gives solitary waves which connect the points O and C~ . 
In Figure 5(a) a -^^2''^ kink with a = 7 = is depicted. The analogous trajectory departs from O and 
arrives at C4. as x goes from —00 to 00. 
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0.05 
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Figure 5: Kink form factor (a), energy density (b) and orbits (c) for K2'^ [h) kinks. 

One of the non-null components, (/)2, is again of kink form but the other, 03 is bell-shaped. The 
energy density is localized at two different points, i.e. this kind of kinks are formed by two basic lumps 
and are two-body solitary waves. They are made from one K^^ and one Kf^ kink, although due to the 
nonlinearity of (|19j) the composition is much more intricate than the simple sum of expressions (|15|) and 
((T7|) . The value of h determines the distance between these lumps; if 6 = 0, the centers of the one-body 
solitary waves coincide |14j and they are superposed. Moreover, 



E[K^^{b)] = E{KY^) + E{Kf^) = 




The energy depends on the eccentricity of the ellipse 63. We remark that the fact of finding the combi- 
nation of one K^^ and one K^'-" as static solutions of © means that there are no forces between these 
basic particles when they are at rest or their relative velocity is zero. The singular member X2''^(0) is 
located on the 02 axis between the points C+ and C_. 
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2. K^^{b): On the ellipsoid E = (Pj + ^ + ^ 
Hamilton-Jacobi equation reads: 



n = ; ^ h 



1, or A = 0, the Hamiltonian H14[) and the 

■ (20) 



dS_ ^ fdS_ dS_ 







Separation of variables S = Sx{x) +S\{X) affords the solution of (|2fl)) for the Hamilton principal 

function via quadratures: Sx = —Ex, 



{F + \Eii)iJi 



^ ■ S 



di>\ 



{F+\Ey)v 



(1 - ii){gI- ii){al- [i] 

if E and F are the separation constants. The finite action trajectories, E = F = 0, are determined by: 
1. Orbits 



dS 



72 



dF 

2. Time schedules 



44> 



1 

'71 



71 



dS_ 

dE 



(1-^)' 



1 



u — ao 



1 



1 



In the field theoretical problem ,the time schedule becomes the kink form factor. 
It is convenient to rearrange the above expressions and write: 



' 1 






" 1 - Z/ " 



















,V2a^(a^y2-2x) 



' 1 


- ^J■ ' 




" 1 - Z/ " 




-ail 















(21) 



Choosing s^ = s^, in ()21() . solitary waves living in in the topological sector C*^*^ are obtained. Back in 
Cartesian coordinates we find: 



1 + g2v^o-2^ _|_ ^ 52 g2v^o-2 

6 



-''32 



52 _|_ £2 g2\/2CTf23; _|_ £^ g-2V2o-|x 



(22) 



1 ^ g-2v^cr|a; _|_ f|. 52 g-2V2(T|2^ 



a, /3 = 0, 1 determine the quadrant of the ellipsoid delimited by the 0i = and (p^ = planes where the 
solutions are confined, see Figure 6(c). Again, this is a one-parametric family of solitary waves, denoted 
as K^^{b). An a = /? = member of this family is depicted in Figure 6(c), asymptotically connecting 
the points Z)+ and -B+. Note that these solitary waves have two kink-shaped non-null components and 
the third one is bell-shaped. The energy density, however, shows two basic lumps, one K^'-" and one 
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Figure 6: Form factor (a), energy density (b) and orbits (c) for -ftT^^ kinks. 



Ki , a composition pointing to lack of interaction between these particles at rest. The energy of these 
two-body kinks depends on the eccentricity of the ellipse 62, 

= EiK^^) + E{K{ 



E[Ki^( 



CD^ _ f^i(2-^^i) 



2V2 



The K^^{0) kink orbit is the elhpse 62; for 6 = E{Kf^) and E{K^^) are superposed, see [Tl] . 

In sum, there are two kinds of static two-body solitary waves formed by either a K^^ plus a 
kinks or a K^^ plus K'^^ lumps. We remark that there is no static configurations involving a K'y^ and 
K^^ kinks or two particles of the same kind. Thus, these particles interact one with each other. 

3.3 Three-body solitary waves 

1. K^^ {a,h):T\ie Hamilton-Jacobi equation 



dS_ ^ fdS_ dS_ dS_ ^ 

dx \dy dn^ dv'' ' ^ 







for the complete Hamiltonian (|14j) of the analogous mechanical system is in any case separable. Writing 
the Hamilton principal function in the form S = —Ex + S\ + + S^, one is led to the quadratures: 



5a 






+ F1A + F3 




A2 

+ T 


2(1 


- X)iai - X)ial - 


A) 




En^ + Fin + F2 






2(1 


-H){aj-fi){al - 






Eu'^ + Fiu + F2 




1/2 


2(1 


- ^){^2 - ^)(^3 - 





Again, finite action trajectories E = F\= F2 = ^ are determined by: 
1. Orbits, given by the intersection of two surfaces in -P(O). 
la. 72 = or. 



(1- A)^(a|-A)4 



(^I-m)" 



(1 - [iy^{[i - a\Yl 



(1 - vY^i{v - a\Yl 
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lb. 73 = H, or, 



(1 - A) ''2<'3 - A) "3^5 



°-32 1 

(1 - - al)'^3 



V^2'l (v — cr|) "2*2 



2. Time schedules, or kink form factor: 71 = or 



_(1-A)^(ct|-A)^ 



tTQ2 3 



These expressions are conveniently rearranged in the simpler form: 

p2V2<72(x-(2-(7|)72+<?i73) 



" A ■ 








V 




1- A 








1 - V 





A 






Sm 


V 


s^ 


Ui-Aj 












" A ■ 


Sa 




Sm 


V 


s^ 















(23) 



„2V2(72(x-(2-(72)72+a273) 



Choosing sa = s^ 



family of solitary waves, denoted as K^^{a,b): 
4>i{x;a,b) 



1 in (|23() we obtain, back in Cartesian coordinates, the two-parametric 



1 + llll e^V2x + f| 52g2^a|x ^ ll a2g2^^2x 
'^3 



^2 



(/)2(x;a,6) 
(^3(a;;a,6) 



^ + g2V2aP ^ 1 g-2V2,7|x ^ 4! 52g2V2,7|2X 

<72 CT^^I ^2 ^2^2 



(24) 



^ + ^ g2V2ais + 1 g-2V2a|x ^ ^1 „2g-2V2ai2X 



Without counting the center of the kinks, related to 71, these kinks depend generically on two param- 
eters a and 6, defined respectively in terms of 72 and 73 as: o"32a^ = e~^^'^2°"32 72^ '7326^ = e~^^'^3°'3273 ^ 
Additionally, a = or 1 specifies in which semi-space 0i > or 0i < these solutions live, see Figure 
7(c). In Figure 7(a) a kink orbit is plotted for a = 0. The particle leaves the point -D+ at x = —00 
and arrives in O at x = -|-cxd, tracing a orbit confined to be in a octant of the interior of the ellipsoid 
E. The associated solitary wave thus lives in the C^^+ sector of the configuration space and has one 
kink-like and two bell-shaped components. The energy density of a solitary wave like this is composed 
of three lumps: the basic particles K^^ ^ K^'-" and K^^ , see Figure 7(b). 

The energy of any kink solution of the 1)24^ family is topological and does not depend neither on the 
(a, b) parameters nor on the coupling constants: 

E[K^^{a, b)] = \W{(i)i{oo; a, b), ^2(00; a, b), ^3(00; a, b)) - W{(j)i{-oo; a, b), (j)2{-oo; a, b), (psi-oo; a, b))\ 

= E{Kn + mn + EiK^) = ^ 



shows degeneracy in energy and a curious kink energy sum rule. 
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Figure 7: Kink form factor (a), energy density (h), and orbits (c) for {a,b) kinks. 

The choice of the parameters a and b to draw Figure 7 aims at identifying the basic components of 
any K^^{a,b) kink. In general, the parameter a is set in order to find kinks in the 03 = plane for 
a = 0, whereas the parameter b runs through M. In Figure 8(c) several orbits of the K^^{0,b) kink 
sub-family are shown. Apparently these solutions are made of two lumps, as illustrated in Figure 8(b). 
The paradox is ficticious because the K^'-" and K^^ lumps are superposed for all the members of this 
sub-family. Therefore, these kinks also involve the three basic particles and we recognize a as the measure 
of the distance between K^'-' and K^^ . 




.25 
0.2 

0.15- 
0.1- 

05 ■ 




Figure 8: Form factor (a), energy density (b) and orbits (c) for {0,b) kinks. 

In the same vein, K^^{a,0) solitary waves live in the plane 02 = 0. Each member of the 6 = family 
displays two lumps, one K^^ kink and a superposition of a ETj^'-" and a K^^ basic particles, b thus 
measures the distance between the ET^*^ and K^^ basic kinks. There exists a member of the family in the 
intersection of the two sub-families, i.e., the K^^[0,0) kink, whose orbit lies on the (pi axis. Generically, 
the K'i{a,b) solitary waves are three-body kinks, but in the 7^3(0,0) configuration the three particles or 
lumps are completely fused. In Figure 9 the energy density of -fC3(10, 6) solitary waves is depicted for 
several values of the parameter b. A similar pattern would be observed by setting the value of b to be 
constant and letting the value of o vary, although the role of the basic particles would be swapped. 



-30 -20 -IC 



-30 -20 -10 



-3 -2 -10 



0.1- 



Figure 9: Energy Density for K^^{1Q, b) kinks: a) b ~ 0, Green zone B, b) b — 5, Green zone B, c) b — 10, Blue 
zone D, d) b ^ 240, Green zone G, and e) b — 10000, Green zone C. All the colored zones are shown in Figure 10. 



It is interesting to remark that the choice of signs, providing the solutions (jSH), is such that the 
inverse image of 

s\Sx{\) + s^S^iji) + SySuif) 
in the coordinate transformation (^HJ is the Hamilton characteristic function in M^: 

1 

71 



1^(01, 02, 03) = ^ \{(tA + 0i + 0^ - 1)' + 2ai02 + 2ai0i] 



(25) 
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Thus, 



1. are the gradient flow hnes of (j)2, ^s): 



dW 



dx 



dx 



dW 



dx 



=F- 



dW 



(26) 



The ODE system (|'26() is separable and tantamount to l\2'6\i in elhptic coordinates. 



2. W{(pi,(f)2,4'3) is no more than the potential energy of a very well known integrable system: the 
Garnier system, the analogous mechanical system of the field theoretical model discussed in j23j 
and [24j . Our model is chosen in such a way that the potential energy of the mechanical analogous 
system is the square of the norm of the gradient of the potential energy of the Garnier system. No 
wonder integrability is found!. 

In sum, there exist three-body solitary wave solutions that are static configurations formed by the 
three kinds of basic particles, with no forces between lumps. 



3.3.1 The moduli space of three-body kinks 



B 



C --10 



I 1- 

-15 -10 -5 



5^^ ^^^^ ^ 

' I:: ^ 



Figure 10: Tiling of the moduli space according to the number of lumps. 

To make these statements more precise, we need a better understanding of the physical meaning of the 
(a, b) coordinates parametrizing the moduli space of solitary waves. To achieve this goal, the number 
of lumps should be identified with the number of maxima of the energy density of the solution on the 
line. For instance, in the (t| = 2ct| = | case, the energy density reads: 



£k{z) = 



2z{2a^ + (3a'' + 4b'^)z + lla^b^z^ + (lOa^ + a^b^ + 12b'^)z^ + (326^ + 2a^b-^)z^ + (27 + &a^b'^)z^ + Sa'^z'^ + b^z'^) 



27{l + a^z + b^z^ + z^ 



where z = e a . It is not possible to find analytically the extremal points of £k{z) because one should 
find the roots of a sixth order polynomial as a function of a and b. Solutions of the equation 



dz 



{z, a,b) = 



(27) 



however, can be studied numerically. We find the following pattern in different ranges of a and b: 



There is a single real solution of ()27|) in the A domain of the moduli space, see Figure 10. The 
solution is a maximum of £Kiz, a, b) and the three lumps are glued together. 
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• In the two B domains there are three real solutions to equation (|27)) ; two of them are maxima and 
the remaining one is a minimum of £k{z, a, b). All over these regions of the moduli space K^"^ and 

sit on top of each other, whereas Kj^^ is apart. 

• Also in the two C domains there are three real solutions, two maxima and one minimum, but now 
only K^^ is separated. 

• Finally, in the four D domains there are five real solutions of (|27|). three of which are maxima and 
the other two are minima of the energy density, and the three basic kinks are totally split apart 
from each other . 



3.4 Four-body solitary waves 

1. K^^{b): The sign combinations s\ 7^ s^ in the (|18j) system of equations offer a one-parametric family 



of solutions belonging to the C^^ topological sectors and confined to living in the 
Cartesian field profiles are: 



plane. The 



(j)i{x) 
h{x) 







-l)°a2a32(l + eV 



(^2 + 0-|2e''e"2y2<7|5; _^ ^^gOg 



32 



a. 



2 b^-2V2alx ^ ^2g-2y2a|xa 



0-3 0-32(1 - e 



-2y/2alx\ 



Again, 6 G M is a real parameter and a = 0, 1. In Figure 11(a) a Kf^{b) kink is represented: the second 
component is bell-shaped but (^3(x) is formed by two kinks. The kink trajectory departs from the point 
B- and arrives in when x goes from — cxd to cxd, crossing through the focus Fi of the ellipse ei. Indeed, 
this behavior is general and the trajectory of each member Kf^{h) passes through Fi at the same "time " 
a; = 0, see Figure 11(c). In Figure 11(b), the energy density is plotted, showing the structure of this kind 
of solitary wave solutions: these kinks are composed of three lumps; two of them are easily recognizable 
as the basic particles K'^^ and Kf^". The other kink is located between the previous ones and has the 
same energy as the sum of the energy of the two basic particles. Later stability analysis will show the 
lump in the middle decaying to the ^ and K^'-' one-body kinks, but exactly superposed in {Q) 
they are in unstable equilibrium. In any case, the K^^{b) solutions are made of four one-body kinks: 
two K^^'s and two K^'^'s. The energy of these solutions is: 



^[Kf ^(6)] = E{KfO) + E[K^^{0)] + E{Kf^) 
An interesting point should be noticed: we claimed that the presence of two particles of the same kind 
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Figure 11: Kink form factor (a), energy density (h) and orbits (c) for K^^ kinks. 
necessarily involves a force between them. Here, we are dealing with a static configuration of two pairs of 
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basic particles. Thus, a delicate balance between the K^^-K^^ and Kf^-Kf^ forces must be reached. 
This is only possible if one K^^ and one K^^ kink share exactly their centers in between the other 
two basic particles, see Figure 11(b). Any small perturbation of this situation would release forces, an 
argument for a dynamical explanation of instability. The fate of these solutions as time passes is explained 
in sub-Section 4.2 of Reference |14j for similar quadruple solitary waves arising in that model. 

2. K'^^{h): Another one-parametric family of kinks exists on the ellipsoid E (A = 0). The sign 
combinations s^ 7^ Sy in the system of equations (|21|) provide the solutions in the original fields: 



4>i{x) 



(-l)"'T2(J3(l + e^e2v^'^i 



{a. 



2 + f^se^e 



32« + o-fge^v^'^i 

0'20'20'32' 



32-^ -I- a^2^°e 



2 „b„~2V2a^x\ 



-1 + e- 



-2V2aix 



(cj2 + (Tie 



-2V2cr^i 



-2sf2a%x 



-l)'3a3a3a32(l + eV2v^-i^ 



(f^32 + ^2^ 



2„fe„-2v^a2 



2„bp2V2a^, 



-|- 1732 e°e 



bp-2V2a^x) 



Two components are bell-shaped, 0i and ^3, whereas the third one, (f)2, is of two-kink form, see Figure 
12(a). b = —^[20^0^212 is the parameter that characterizes the family members; the values a, /? = 0, 1 
specify the quadrant of the ellipsoid E in which the kink is located. These solitary waves belong to the 
C'^'-' topological sectors and connect the points C"*" and C~ crossing the umbilicus points see Figure 
12(c). In Figure 12(a), a K^^ [b) kink is depicted leaving from point C_ and arriving at C-|_ as x goes 
from —00 to 00. In Figure 12(b) the energy density is shown, displayed by three lumps; two of them are 
easily identified as the basic particles K^^ and Kf^^ but the more energetic lump, located between the 
previous ones, corresponds to a combination forming the K2^{Q) kink. Thus, the K^'-'{b) solutions are 
composites of four basic particles, two K^^ and two K^'-' ones. The energy depends only the coupling 
constant CJ3: 



al{2 



V2 



The K^^ {b) four-body kinks are qualitatively identical to the Kf^{b) kinks replacing the ellipsoid E 
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Figure 12: Kink form factor (a), energy density (h) and orbits (c) for Kf^ kinks. 

by the 4>i = plane (A = and v = 1 faces in the P{0) cuboid). Thus, all considerations about stability 
and forces made for the Kf^{b) kinks also work for the K^^ {bys. 



3.5 Six-body solitary waves 

1. K^^ {b): There is a more complicated one-parametric family of solitary waves in the C^^ topological 
sectors living in the (;/>3 = plane. The potential is restricted 

02,0) = (</.? +0^)(0?+</.i 



lf + 2al(t>l{ 



^1 + ' 



1) + at4>l 
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to that discussed in the two scalar field model of Reference Jl]. In the present context, however, the 
Hamiltonian ()14() reduces to either the I) A 



(t| or II) /Lt = cj| faces of P(0): 



Hi 



+ 



(A 



+ 



The Hamilton-Jacobi equations for Hi and Hn must be solved consecutively and the two pieces should be 
continuously glued afterwards, see Reference [53] for similar solutions in a simpler model with a quartic 
potential energy. Sign combinations for sa, s^, s,y leading to these solitary waves are different from those 
giving the 7^^3*^(0,6) kinks. 

These solutions connect the points and C~ and cross the F2 foci (see Figure 12(c)). The kink 
form factors are: 

{-\Ya2 fl + e^v^-^^ 



al + a2e-f'e-2v^^ + e2v^-i^/l + ale^^'-l^ + ale'' e-'^^'^l'^ 



^2V2alx + ^2^2V2x + ^2 b Lb ^2V2alx ^ ^2^b + ^2g2^/2x 



())l{x) = 

02 (x) = 
<^3(X) = 

In particular, a K^'~" {b) kink with a = has been depicted in Figure 12(a), departing from the point 
C- and arriving in C+ when x goes from —00 to 00. One component is of bell form and the other is 
two kink-shaped. In Figure 12(b) the K^^ {b) energy density is shown. Apparently, these solutions are 
formed by three lumps. One of them is the basic particle K^^] the more energetic lump in the middle 
of the other two is the 7^3*^(0) kink (the union of the three basic particles), and the third lump is a 
K2'~' kink (superposition of the K^^^ and JC^*^ basic kinks). Therefore, the Kq^ {b) solutions involve 
six basic particles, a pair of each type of basic lumps. These static configurations describe an even more 
delicate equilibrium than four-body lumps. The forces between basic kinks of the same type are only 
balanced at the expense of joining them in a arrangement such as that shown in Figure 12(b). Any 
relative displacement of the centers of the composite particles would destroy these unstable solutions. 
The energy of these six-body solitary waves does not depend on the coupling constants: 



E[K^^{b)] 



E{K^^) + E[K^^{{), 0)] + E[K^^ 



2E{K^^) + 2E{Kf^) + 2E{Kf^) 
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Figure 13: Kink form factor (a), energy density (h) and orbits (c) for K^^ kinks. 



2. Kq : Another one-parametric family of six-body solitary waves exists in the C topological 
sectors, now located on the plane (j)2 = 0. Also, the restricted potential energy, 

V{<Pi,0, ^2) = {<lA + 0i)(0? + <fi- if + 2ai0i(0? + + - 1) + atcfi , 
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corresponds to a two-component scalar field model studied in jj^. There is one difference in this three- 
component situation: as in the previous case, the Hamiltonian ()14() reduces to either the I) ^ = ^2 or II) 
i' = a2 faces of P{0): 



Hi 



+ 



H. 



{X-ai){\-v) {u-X){u-ai) 



Hii 



+ 



H„ 



in - A)(/i 



Again, the Hamilton-Jacobi equations for Hi and Hn must be solved consecutively and the two pieces 
must be continuously glued afterwards. The choice of sign combinations for sx, s^, is opposite to that 
giving the /rf^(a,0) kinks. 

These solutions asymptotically join the points and crossing the foci F3. The kink form factors 

are: 



1)°CJ3 (1 + e^v^-^i^ 



(j)l{x) 



In Figure 14(a) a K^^ kink is depicted with a = that connects the point B- with B^. The structure 
of one member of this family is similar to that found for the previous kinks. From the energy density 
plot shown in Figure 14(b), one sees that each member of the family consists of six basic particles: two 
, two K^^ and two K^'-^ kinks. Only the order differs (the two-body kink is on the right of the 
three-body lump) and the two-body kink is made of a K^^ and a K^'-" solitary wave, instead of a ^ 
and a Kf^. The energy, however, is the same: 



£;(i^g««) = ^(Kf ^) + E{K^^{0, 0)) + S(i^2 ) = 2S(i^i^^) + 2E{KY'') + 2E{Kf'-) = ^ 
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Figure 14: Kink form factor (a), energy density (h) and orbits (c) for K^^ kinks. 



3.6 Seven-body solitary waves 

2. Kf^{-i2,l?,): Besides the generic (depending on two parameters) K^^{a,b) kinks, there are other 
generic solitary waves. Other combinations of signs in the system of equations (|2,'-{)) provide new generic 
kinks different from K^^ {a,b). The kink orbits of these solutions run between the B and C minima 
tracing a curve with three well-defined stages: 1) First, the kink trajectory departs from the point C 
at x = — oo remaining in a octant of the ellipsoid E interior until the particle reaches a point of the 
characteristic hyperbola h. During this stage the trajectory satisfies the system of equations (|2^-i|) with 
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sa = 7^ Si,. 2) Then, the kink orbit crosses to the next octant (from left to right or viceversa) through 
this point and traces a curve that ends by hitting a point in the characteristic elHpse 64. In this second 
stage the kink trajectory also complies with (|2l-ij) but the sign combination is s\ = 7^ s^. 3) The last 
stage of the kink trajectory starts again by changing the octant (from up to down or viceversa) through 
the point in 64. The kink orbit runs along a curve in the latter octant until the point B is reached. 
The final stage satisfies the system with sx ^ = s^, see Figure 15(c). The parameters 72 and 
73 respectively determine the points where the kink orbit hits the characteristic hyperbola h and ellipse 
64. Orbits for which the relation 72 — (T373 between integration constants holds meet at a single point 
on ellipse 64. The same happens for trajectories that comply with 72 — ct|73 meeting at a point on 
the hyperbola h. Thus, at each point in h and 64 an infinite number of kink trajectories meet. These 
characteristic curves are focal lines. 

These solutions are so complex that it is not possible to invert the relationships (|23j) in order to 
obtain analytical expressions in Cartesian coordinates for the kink form factors. For this reason, we 
parametrize these solutions in terms of the integration constants 72 and 73 that set the orbits in the 
Hamilton- Jacobi expressions H23() . The use, however, of the implicit expressions in elliptic coordinates 
allows Mathematica to plot field profiles and energy densities. It should be mentioned that extreme care 
is necessary to guarantee continuity of the trajectories and its derivatives when it rebounds on the faces 
of the -P(O) cuboid. In Figure 15(a) the form factor (field profile) of a ^^'^(72,73) kink is shown. In the 
Figure 15(b) we display the energy showing seven basic lumps: two Xf^'s, two K^^^s and three X^'^'s, 
with the peculiarity that three different basic particles always travel together. We remark that there 
exist no net forces between particles in these configurations but that they are highly unstable against 
any perturbation splitting the centers of superposed particles. The total energy is: 



2E{K^^) + 2E{K\"') + SEiKi"^') = ^ + 
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Figure 15: Kink form factor (a), energy density (h) and orbits (c) for kinks 



A summary of all the kink energies is offered in the following table: 



Type of Kink 


Energy 


Basic Kinks 




Two-Body Kinks 




Three-Body Kinks 




Four-Body Kinks 


E[Kf^{b)] = % E\K2^{h)\ = ^-^^ 


Six-Body Kinks 


EiK^'^m = EiK^ib)] = ^ 


Seven-Body Kinks 


E[K-c^a,b)] = i, + '-l^ 
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